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Abstract 

We prove the global well-posedness of three dimensional compressible Navier-Stokes 
equations for some classes of large initial data, which is of large oscillation for the density 
and large energy for the velocity. The structure of the system (especially, the effective 
viscous flux) is fully exploited. 

1 Introduction 

We consider the three dimensional compressible Navier-Stokes equations 
dtp + div(pu) = 0, 

d t (pu) + div(pu <g> u) - pAu - AVdhm + VP(p) = 0, (1.1) 
(p,u)\ t =o = (po,u )- 

Here p(t, x) and u(t, x) denote the density and velocity of the fluid respectively. The Lame 
coefficients p and A satisfy 

p>0 and A + /i>0, (1.2) 

which ensure that the operator p Au + A Vdivu is elliptic. The pressure P is a smooth function 
of p and satisfies 

P(0) = 0, P'(p) > 0, (1.3) 
(p-p)(P(p)-P(p))>0 for pj^p,pe [O^Cq 1 ], (1.4) 

where < cq < ~p < Cq are fixed numbers. 

The local existence and uniqueness of smooth solution of the system (jl.ip was proved 
by Nash [31] for smooth initial data without vacuum. In a seminal paper [30], Matsumura 
and Nishida proved the global existence and uniqueness of smooth solution for smooth initial 
data close to equilibrium. In general, whether smooth solution blows up in finite time is an 
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open problem. However, Xin [34] proved that smooth solution will blow up in finite time if 
the initial density has compact support. Recently, Sun, Wang and Zhang [33] showed that 
smooth solution does not blow up if the concentration of the density does not occur. 

The global existence of weak solution was proved by Hoff [22j |25] for the discontinuous 
initial data with small energy. For the large initial data, the global existence of weak solution 
was proved by Lions [29] for the isentropic Navier-Stokes equation, i.e. P = Ap 1 for 7 > |. 
Feireisl, Novotny and Petzeltova [20] improved Lions's result to 7 > \. For the spherically 
symmetric or axisymmetric initial data, Jiang and Zhang |27[ [28] proved the global existence 
of weak solution for any 7 > 1. However, the question of the regularity and uniqueness of 
weak solution is completely open even in the case of two dimensional space. 

Motivated by Fujita and Kato's result for the incompressible Navier-Stokes equations [21], 
Danchin in a series of papers |15j fT6"j \T7\ [TO] proved the local well-posedness of for the 

. i . 3—1 

initial data (po — ~p,uo) belonging to the critical Besov spaces B^ x B£ 1 for p < 6, and 

.3 . \ 

global existence for the initial data close to equilibrium in S| 1 x B^ \ ■ This seems the largest 

space in which the system (II. ip is well-posed. Indeed, Chen, Miao and Zhang [TJ] proved the 

. 3 . 3_ 1 

ill-posedness of (jl.ip in B'jX B^ x for p > 6. 

Cannone, Meyer and Planchon [3j H] proved the global existence of the solution for the 

incompressible Navier-Stokes equations in the Besov space of negative regularity index Bp oo 
for p > 3, see also pQ for a new proof based on the smoothing effect of the heat equation. A 
very interesting point of this result is that it allows to construct global solution for the highly 
oscillating initial velocity that may have a large norm in or L 3 . A typical example is 

u (x) = sm(—)(-d 2 (i)(x),di4>(x),0) 

where <fi G 5(R 3 ) and e > is small enough. Later, Chemin, Gallagher, and Paicu [9] 
find more classes of initial data with large energy generating the global solution. And in the 
proof of [9], the special structure of the equation is used. We refer to [ITJ [8j [lOj [32] [36] and 
references therein for more relevant results. 

It is highly non trivial to generalize Cannone-Meyer-Planchon's result to (jl.ip . since it 
is a hyperbolic-parabolic coupled system. Very recently, important progress has been made 
by Chen, Miao and Zhang [13] and Charve and Danchin [5], where they construct the global 
solution of (jl.ip for the highly oscillating initial veloicty by proving the global well-posedness 
of the system (jl.ip in the critical Besov space with p > 3. 

However, for all the above global existence results, the initial density is required to be 
close to a positive constant in L°° norm, which precludes the oscillation of the density at 
any point. Recently, Zhang [35] proved the global existence and uniqueness of (jl.ip for the 
initial density po close to a positive constant in 1? norm and uq small in LP norm for p > 3, 
which allows the density to have large oscillation on a set of small measure. Similar result 
has also been obtained by Huang, Li and Xin [26] for the initial data with vacuum, but a 
compatibility condition is imposed on the initial data. 

The proof of [35] and |26j is based on Hoff's energy method, where the structure of the 
system is fully exploited, especially the effective viscous flux plays an important role. While, 
the proof of |13] and [5] is based on the analysis for the linearized system, where Harmonic 
analysis tools (especially Littlewood-Paley theory) are exploited, but the structure of the 
system is not exploited. 
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In this paper, we combine two methods to find a new class of large initial data generating 
the global solution of the system (jl.ip . This class of initial data may have large oscillation 
for the density and large energy for the velocity. More precisely, we prove 

Theorem 1.1 Assume that the initial data (pq,Uq) satisfies 

Po-pGiTni^, c < Po <c \ n G H s - 1 n B*' 1 (1.5) 
for some p G (3,6) and s > 3. There exist two constants c\ = ci(A, pi, Co, p) and 

c 2 = — —~ m — C = C ( X >f J '> c 0'P) 

C(l+ p -p 3 + 1/0 2 



,i 



snc/i t/tat i/ 



P0-P||l 2 < c 1> Ikollff-S + 1 1 T*0 1 1 3-1 < C2, (1-6) 



p 

p,l 



/or some <5 G (1 — |, |), t/ien t/iere exits a unique global solution (p,u) to the compressible 
Navier-Stokes system hl.l\) satisfying 

P>f, p-peC([0,+oo);H s ), 
u G C([0, +oo); £P _1 ) n L 2 (0, T; fP) /or any T > 0. 

Remark 1.2 Since the constant c\ is independent o/po, i/ie condition ||po — p||l 2 — c i allows 
the initial density to have large oscillation on the set with small measure. Given po, if we 
take the the highly oscillating initial velocity, for example, 

x 

uq(x) = sin(— )<p(x) 
for any ip G 5(R 3 ), then for p > 3 and e small enough, 

,ii ii ^ n min(<5,l — -) - 

+ | «o|| . 3-i < Ce y p' < 



UO i-S + PO 3_! S 0£ p S c 2 . 

OP 

c p,l 



Interestingly, uq has a large norm in H s for s > 0, hence the large energy. 

Remark 1.3 Whether similar result can be generalized to the case of two dimensional space 
remains unknown. However, it seems impossible to generalize this result to the full Navier- 
Stokes equations, see [Lffl . 

The main steps of the proof are as follows. First of all, we construct a local solution of (jl.ip 

. 3 . 3_ 1 

for the initial data (po — p, no) belonging to the critical Besov space B^ x x B* 1 . To deal with 
the case of physical non vacuum (i.e., po > Co > 0), we need to use the weighted Besov space 
introduced by Chen, Miao and Zhang [12]. Using the smoothing effect of the momentum 
equation, this solution can propagate the smallness of the initial velocity in Sobolev space 
of negative regularity index to that of the velocity at some time to > in Sobolev space 
of positive regularity index. Then Hoff 's energy method can be applied to prove that the 
density is bounded below and above after t > to- Finally, we prove a continuation criterion of 
smooth solution, which ensures that the local solution can be extended to a global solution. 
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(2.1) 



2 Littlewood-Paley analysis 

2.1 Littlewood-Paley decomposition 

Let us introduce the Littlewood-Paley decomposition. Choose a radial function ip G 5(R 3 ) 
supported in C = {£ G R 3 , § < |£| < §} such that 

^^(2-^) = l for all £^0. 
jez 

The frequency localization operator Aj and Sj are defined by 

A j f = <p(2-W)f, Sjf= ^ A k f for jGZ. 

k<j-i 

With our choice of ip, one can easily verify that 

AjA k f = if \j - k\ > 2 and 
A j (S k - 1 fA k f) = if U-fc|>5. 

Next we recall Bony's decomposition from [2]: 

uv = T u v + T v u + R(u,v), (2.2) 

with 

T u u = 52 Sj-iuAjV , iZ(ti, v ) = Aj«AjV, Ajv = 52 ^j' v - 
jez jez \j'-j\<i 

The following Bernstein lemma will be frequently used(see pQ). 

Lemma 2.1 Let 1 < p < q < +00. Assume that f G L P (R 3 ), i/ien /or any 7 G (NU{0}) 3 , 
i/iere exist constants C\, C2 independent of f, j such that 

supp/ c {|e| < 4,2*'} < C^+W^WfWp, 

supp/ C < |e| < A 2 2>'} ll/llp < C 2 2-M sup ||^/|| p . 

|/3|=M 

2.2 The functional spaces 

We denote the space Z'(R 3 ) by the dual space of Z(R 3 ) = {/ G <S(R 3 ); D a f(0) = 0; Va G 
N 3 multi- index} , which can be identified by the quotient space of S'(R 3 )/-P with the poly- 
nomials space V. 

Definition 2.2 Let s G R, 1 < p, r < +00. The homogeneous Besov space B^ r is defined by 

b; = {/^'(r 3 ): ii/Hb* <+oo}, 

where 

11/11^ ^ 2 fcs ||A fc /(t)|| p . 
Ln particular, -Bf 2 = H s , H s is the homogeneous Sobolev space. 
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Let {e k (t)} k€ z be a sequence denned in [0, +00) satisfying the following conditions 

e fe (t) € [0, 1], e k (t)<e k >(t) if k < k' and e k (t) ~ e k >(t) if ~ k', (2.3) 

where ~ k' means that there exists a constant iVo such that \k — k'\ < Nq. The weight 
function {co k (t)} k& z is defined by 



U) k 



It is easy to verify that for any k £ Z, 

Wfe(i) < 2, e fc (i) < Wfc(t), 

^fe(i) < 2 k ~ k 'uj kl (t) if > fc', < 3o; A ./(t) if fc < k', 

<^>k(t) ~ if ~ k! . 

The following weighted Besov space is introduced in |12j . 



(2.4) 



Definition 2.3 Le£ s G R, 1 < p, r < +00, < T < +00. The weighted Besov space Bp r (u) 
is defined by 



where 



d = f \\2 k °uj k {T)\\A k f\\ p \\ tT . 



Obviously, Bp r C Bp T {uj) and 



Bi,r(w) 



< 2 



We also need to use Chemin-Lerner type Besov spaces introduced in [6]. 

Definition 2.4 Let s£R, 1 < p, q,r < +00, < T < +00. The functional space Lj,(Bp 
is defined as the set of all the distributions f satisfying 



p,r/ 



def 



P,r) 



2 fes ||A fe /(t)|| 



< +OO. 



The weighted functional space L^p(Bp x(u;)) is defined similarly, whose norm is given by 



def 



2 ks u k (T)\\A k f(t)\\ 



2.3 Nonlinear estimates in Besov spaces 

Let us recall some nonlinear estimates in weighted Besov spaces from [12 
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1,1 1 

i y> yii ^2 ^ w*" 1 



Lemma 2.5 Lei 1 < p, q, q\, q^ < oo with ^- + ^ = i TTien £/iere ZioZd 



l|Tff/ll z«(B; i 1 + ' a -^( W )) - c ' ll/ll ^ 1 (^ 1 (-)) |l5llz T 2 (s;f 1 ); 

f&j i/sx < | — 1, we have 

l{Tf9 huB p \ +32 -h»)) ~ c|I/IIz t 1 (^ 1 h) i|5|I ^ 2 (^ 2 1 ) ; 

(cj if s\ + S2 > 3max(0, | — 1), we /iawe 

The following lemma is a direct consequence of Lemma 12.51 

Lemma 2.6 Let si < | — 1, S2 < |, si + «2 > 3max(0, | — 1), and 1 < p, (/, gi, g2 < oo wrai/i 

— + — = - . Then there holds 
qi 92 g 

ll/^IL, (B ^ 2 -I H) ^ cil/llz^^c-)) 11 ^ 1 ^ 2 ^)- 

In general Besov spaces, we have 

Lemma 2.7 Let s±, S2 < |, si + S2 > 3max(0, | — 1), and 1 < p, qi, q-z < oo u>i£/z — + — = 
-. T/ien i/iere /io/ds 

9 



H/5II n+S2 -| < C\\f\\^ T HB;iJ9\\l ¥iBZ y 



L q T {B p \ p ) 

Lemma 2.8 Let s > and 1 < p,q,r < oo. Assume that F E W^ +3 '°°(R) with F(0) = 0. 
Then there holds 

\\ F (f)huB^)) * C{1 + ll/ll^(^)) H+2 H/llz^(^H)- 

T/ie same result also holds true in Besov space without weight. 

The following is a weighted commutator estimate. 
Lemma 2.9 Letp^[l,oo) and s G (— 3min(|, ^-), |]. Then there holds 

11^(2011^,^11^,1^ < CH/II a l| 5 |lzi ( ^ V) . 

L r ( B p ,iM) p ' 

Proof. Using Bony's decomposition (|2.2p . we write 

[/, Aj] • <9 fc5 = [T f , Aj]d k g + T dkA . g f + R(f, d k A j9 ) 



Using Lemma [231 (a) and (c) with s\ = | and S2 = s, we get 



H^'WAiiTa^m^L <cil/ll ~ f j , Jak^r 

|^(T)2^||A J ii(/ fe ,5 fc5 )|| i i (LP) || £1 < C||/|L 3 WgWzufr+iy 



Noticing that 



Ta k A j9 f = Ta k A j9 f + R(f,d k A jg )= ^ S' i , +2 A,a fc5 A / /, (2.5) 

i'>j"-2 

then we get by Lemma 12.11 and (|2.4p that 

||^(W-||I^ Aja /|| r x, (£p) ||^ 

<C||o; i (r)2^||A J -V 5 ||^ (Loc) ^ ||A f /||^ (iP) ||^ 

i'>j'-2 

<C||2 j(s+1+ f)||A^|| LMLP) ^ M^i'/IU^)!!* 



< c||/IL .| iiffiip.fBH-i,- 



j'>i-2 



Set = (J 7 Then we have 

\T f ,A 3 \d k g = ^ / J 'yVS j ,- 1 f(x-Ty)dTd k h@y)A j ,g(x-y)dy 

+2 3 ^' / h(2^x-y))d k S j ^ 1 f(y)A f g(y)dy, 
Jn 3 



from which and Young's inequality, we infer that 

mCT^lCT/, A^lli^lU <C||/|| a ll^llzi^t 1 )- 

The proof is finished by summing up all the above estimates. ■ 

3 The linear transport equation and momentum equation 

We consider the linear transport equation 

d t f + v-Vf = g, 
f(0,x) = f . 

Set V(t) ^ J* * ||Vi>(t)|| 3 dr. The following result is from 

B p,i 



7 



Proposition 3.1 Let p G [1, +00] and s G (— 3min(i, ^7), 1 + |). Let v be a vector field such 

3 

£/ia£ Vt> G L\,{B^^). Assume that /o G -B*i, 5 G L^i?^) and / is i/ie solution of 113.1]) . 
Then there holds for t G [0, T], 



The following weighted version is from [12] . 



Cnt) (\\f \\„s + [ e- cv ^\\a(r)\\ as dr). 



Proposition 3.2 Let p G [l,+oo] and s G (— 3min(-, — ), -]. Let v be a vector field such 

that Vv G L\,{B^^). Assume that /o G g G L^B^) and f is the solution of 113.1]) . 

Then there holds for t G [0,T], 



^ OT(,, (ii/oiiB. lM +r^ w ii 5 wiiB. lM * 



We next study the linear momentum equations with variable coefficients 

dtu — div(/ZVti) — V(Adivn) = G, 
u\ t =o = u . 



(3.2) 



Assume that the viscosity coefficients X(p) and p(p) depend smoothly on the function p and 
there exists a positive constant C3 such that 

7^>c 3 , A + 7Z>c 3 . 

And the weighted function u>k{t) is taken as follows 

u k (t) = Y,2 k - e ee(t), 

e>k 

with ei(t) = (1 — e _c22£ *)2 for some positive constant c > 0. It is easy to verify that the 
function e^(t) satisfies (|2.3p . 

Proposition 3.3 Let p G [2, 00), s G (1 - ~,~],<? G [l,oo]. Assume that u G B s p ~^,G G 
L} r (Bp'] 1 ), and p — ~p G Lj?(B£i)- Let u be a solution of \3.2\) . Then we have 

H n llz^(^+ 1 ) + II^IIZKs^) 

< c(||«o||^- 1(w)) + ll^lll^^H) +^( r )IIP-pll £?( ^| iH) ll«llzM< 1 1 ))' 
and t/ie version without weight 



C lko||^-i + +A(T)\\p- p\\_ 3 Hulljw^+is 

V p,l ^TV^p.l J ra/nP \ ^TV J>,1 ' 



HereA(T)t f (l + \\ p \\ {Lao) )^ +2 . 
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Proof. We only prove the first inequality, the proof of the second inequality is similar. Set 
d = diva and w = curkt. By (|3.2fl . (d,w) satisfies 

d t d - div(z7Vd) = divG + F u 

d t w - div(jlVw) = curlG + F 2 , (3.3) 
(d,w)\ t =o = (divu , curln ) = (d ,w ), 



where v = A + fx and 



Ft = div(V/Z • Vn) + div(V(A + Ji)d), 

F 2 iJ = div^/ZVu* - diJlVu j ), i, j = 1, 2, 3. 



Apply the operator Aj to (|3.3p to obtain 

f <9 t A,,d - div(FVAjd) = divAjG + AjF x + div[Aj,T>]Vd, 
\ d t AjW - dfrijiVAjw) = curlA^G + AjF 2 + div[Aj,Ji]Vw. 

Multiplying the first equation by | Ajd\ p ~ 2 Ajd, we get by integrating over R 3 that 
- — \\Ajdl\P- / diviWAjd^Ajd^Ajddx 

(divAj-G + AjFx + div[Aj,i/]Vd)|Ajd| p ~ 2 Aj(Ws. 



'R 3 

There exists c p > depending on c%,p such that (see [16]) 



This gives that 

d 
~dt 

which implies that 

r-t 



WAjdWp + ^WAjdWp < C(2^||A i G|| p + ||A i F 1 || p + 2^||[A i ,z7]Vd|| p ), 



W&AQWp < e- Cp22Jt ||A^ || p + C f e-^^i^WAjGWp + ||A,iq|| p + 2^'||[A i ,I7]Vd|| p )dT. 

J o 

Similarly, we can show that 

\\A jW (t)\\ p < e-^ 22Jt ||A^ || P + G [\- c ? 223 ( t - T \2i\\A j G\\ p + \\A j F 2 \\ p + 2i\\lA J ,Tl}Vw\\ p )dT. 

Jo 

Hence, we infer that for any q £ [l,oo] and t G [0, T], 

\\&jd{t)h«(LP) + \\ A jw{t)\\ L i tiLP) 
<C2- 2 ^c j (T)k\\A j d \\ p +\\A j w \\ P ) 
+C2-V/%(T) 1 *[2l\\A j G\\ Llm + IIA^H^^ + 2J||[A J -,F]Vd|| i i (iP) ) 

+G2- 2 ^c i (r)i(2^||A i G|| il(LP) + HA^IU^) + tf\\[A j ,]i]Vw\\ LULP) ), (3.4) 
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with Cj(T) = 1 - e- c f 22JT . Notice that 

2 j \\A jU \\ p ~ ||Ajd|| p + \\Ajw\\ p , ej{T) < uij(T), 
which along with (|3.4p implies that 

< Cdluoll^-i^j + l|G|| z ^ ri i H) + IK^i, ^2)|lzi ( ^- a(w)) ) 

+C\\2^w J (T){\\[A j ,u]Vd\\ LULP) + HfA^Vu;!!^))^. (3.5) 

From Lemma 12.61 and Lemma 12.81 it follows that 



H^Hz^-H) + H^Hr^co,)) ^ a(t)IIp_pII ^ ?( bI i(w)) i|u|I ^(^ 1 )' (3 - 6) 

and using Lemma 12.91 we get 



< A( T )\\P ~ P\\~ .1 Nip (B s + l V ( 3 - 7 ) 
Summing up (|3.5p - (|3.7|) yields the first inequality of Proposition! 



4 Local well-posedness with physical non vacuum 

In this section, we prove the local well-posedness of compressible Navier-Stokes equations 
in the critical Besov space with physical non vacuum assumption. That is, we just impose 

_ 3 

Po > cqO for the initial density. In the case when p — ~p is small in BE i or has more regularity, 
the corresponding local- well posedness has been obtained by Danchin [16j. Recently, Chen- 
Miao-Zhang [12] and Danchin [19] developed two different methods to deal with physical 
non vacuum case. Here we will revisit the proof of [12 . The main goal is to record the 
precise dependence of energy bounds and the existence time T on the initial data and other 
important constants, which is very important for our argument. 

Theorem 4.1 Let ~p > and cq > 0. Assume that the initial data (po,uo) satisfies 

3 3 i 

po - p G B^, c <p < c \ u G B* A . (4.1) 

Then there exists a positive time T such that 

(a) Existence: If p G [2,6), the system has a solution (p — p,u) G Ej, with 

i^C^T];^)* (caO.Tj^l^nLHo.r;^ 1 )) 3 , p > \c Q] 



(b) Uniqueness: If p G (1,3], then the uniqueness holds in Ej,. 
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Proof. We only present the uniform energy estimates for smooth solution of , since 

the existence can be deduced by constructing approximate solution sequence and compact 
argument. The uniqueness has been proved in [12]. Set 

a{t,x) = , f Jl (p) = -, A(p) = -. 

P P P 

Then the system (jl.ip can be rewritten as 

dta + u ■ Va = F, 

d t u - div(/ZVu) - V(Adiv u) = G, (4.2) 
(a,u)\ t =o = (ao,u ), 

with a = and 
u p 



F(= -(1 + a)divu, 

G = —u-Vu — — Va + 4- v P • Vu + ^Vpdivu. 

p p z p z 

Let us assume that the solution (p, u) satisfies 

j<p(t,x)<2c Q 1 on [0,T]; (4.3) 
||o|L a + lldl 2_i < Co(||ao|| a + 1 1 i*o 1 1 2-i) = f CoE ; (4.4) 

IkL 3 <vi; (4.5) 

llull 3 + llttll 3 < m- (4.6) 

These assumptions are satisfied for the solution in Ej, if T is small enough. We will show 
that there exists a T > depending on the initial data and A, p,,p, cq,p such that these 
assumptions are satisfied. 

In what follows, we denote C by a constant depending only on A, fj,,~p, co,p, which may be 
different from line to line. 

Step 1. Estimates in the Besov space 

Applying Proposition 13. 1 1 to the first equation of (14. 2p . we get 

[loll 3 < e cy(T) (||ao|| 3 3 ), (4.7) 

and applying Proposition 13.31 to the second equation of (|4.2p . we have 

NL 3-1 < C||«o|| a_ 1 + CIIGII s_! +C|lo[| 3 [lull 3 +1 , (4.8) 

where = J ||n(r)|| 3 +1 dr. By Lemma 12.71 we have 

|IjP]| 3 < [lull 3 +1 + C][o[l a ||it[L i +1 , 
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and by Lemma 12.71 and Lemma [X 

||G|| 3 ! <C\\U\\_ .1.JML 3 +1 +C[|0|U 3 (T+HL 

Here we used p < 6. Plugging them into (|4.7p - (|4.8p . we obtain 

||o||_ a +HI a-i < CeC??2 (^0 + ^0(^2 + T) + r/ 2 ). 

Step 2. Estimates in the weighted Besov space 

Applying Proposition 13.21 to the first equation of (|4.2p . we get 

IkIL .a <e cy ( T )(||oo||.3 +||F|L .3 ), 

and by Proposition 13.31 we have 

ll n ll 1+1 + ll n ll 2 
L^BfJ ) 

<cf||«o|| 3-1 +I|G|L a_i +lk|L 3 \\u\\_ 3 +1 ) 

V Bfr («)) ^(B^ («)) L^(B^)) L^B^ y 

By Lemma 12.71 we have 

II^IL 3 <2||F|L .3 <C(l+||a||_ .2 )\\u\\_ .3 +1 , 
4(^ p i) ^(^1) l^b^ ) 

and by Lemma 12.61 and Lemma 12.81 

||G|| 3_, <C\\uf 3 + C||a||_ 3 (T+IML 

L^B^ («)) ^(%M) ^(^,1 ) 

Plugging them into (|4.1U|) - (|4.11|) yields that 

||a|L .1 <Ce Cl??2 (||ao|| .3 + (1 + C E ) m ) , 

L^B^fyj)) ^i(«) 

IML .1+1 + NL .2 < CfNoll . 3_! +f?| + 77i(T + 77 2 )). 

Step 3. Estimate of the density 
Let x) be a solution of 

x) = u(t, X(t, x)), X(0, x) = x, 
and we denote by X~ l {t, x) the inverse of X(t, x). Then a(t, x) can be solved as 

a(t,x) =a {X- 1 {t,x)) + f F{T,X{T,X-\t,x)))dT, 

Jo 

this gives 

p(t,x) = Po (X- 1 (t,x))+p f ' F{T,X{r,X-\t,x)))dr. 
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Since we have 

II^IL^L-) < l|V'u|| L i i(L o 0) (l + |MUs?(L°°)) < Cf] 2 , 
from which and (|4.14p . it follows that 

p(t,x) > c -Cri2, p(t,x) <c^ + Cr] 2 . (4.15) 
Step 4. Continuity argument 

Let Co = IOC. We first take 772 and 771 small enough such that 

1 



Cm < 9, e Cr » < 2, Com < 1, m< e , 

o 

Cm<j, 2C(i + c E )m < T j- 



(4.16) 



Then we take T small enough such that 



C T<1, T< V2 , \\ao\\ a Nil a-: < (4-17) 

B* tl ( u ) 4C ( W ) 4C 

With this choice of 771,772 and T, it follows from P~9]h (|CI2]> . (f4T3|) and P~T5|) that 

4 

o l 3 + M a_i < tCo-Eo, 

^(^! ) 5 

M H 3 

ML .3 < -rVl, 

L^CB^H) 4 

3 

ML .3+1 + ML .3 < 7772, 

) l|(b; :1 ) 4 

3 5 
-c < p(t,z) < -c \ 

which implies that the assumptions (|4.3|) - (|4.6|) are satisfied for T defined by (|4.17|) . ■ 

Remark 4.2 In general, the choice of time T depends on Co = Co(A, p,, p, co,p) and the 

_ 3 

behaviour of the initial data (po,uo). If \\uq\\ 3_ 1 < c and oq G Bi\ n H s for s > § and 

B p , 

c = c(\, p,~p,co,p) small enough, then T can be taken as 

T = — C = C(A,/ J ,p,c ,p, S ). (4.18) 

C l + ||ao|| 3 V s +2)/U 5) 

Indeed, due to f^.ifip , we can ta&e 771,772 as follows 

1 1 

*7i = ~> ^2 



C " C(l + ||a || .3 ) 



for some constant C depending only on X, p, p,co,p. And by Lemma \2. 11 for any N £ N, 
2^||A J a || P < C ^W^M* < C2-^l) JV ||a ||^. 

j>N j>N 
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And by {2.1$ , we have 



y~] 2 J, 'po; i (T)||A i ao||p < o;jv(T)||oq|| 3 . 

— B p 

3<JV P.i 



Then j^. is ensured by taking T as {4-18 ). 



5 Propagation of Sobolev regularity 
5.1 Propagation of low regularity 

Recall that in Section 4 we have constructed a solution (u, p) of satisfying 

^7 < p(t, x) < 2CQ 1 , llall 3 + [lull 3_! < C E , 
2 ^ ^(^1) ) 

||a|L 3 < r/i, 3 +1 + 3 < 772, 

with a(t, x) = x) — -Eo = Ikoll a + ll n o|| S-i- 

We will show that this solution allows to propagate the regularity of the initial data in 
Sobolev space with low regularity. This will ensure that one can use Hoff 's theory for small 
energy initial data, since the norm of the initial data is small in Sobolev space with low 
regularity, but it is large in Sobolev space with high regularity. 



(5.1) 



Proposition 5.1 Let p 6 (3,6) and 1— | < 8 < -. Assume that (p,u) is a solution of 

rl-S ^ ii-8 



satisfying \5.1\) . If (a ,u ) € H x H °, then 

W a h™(Bl- s ) ^ C {\\ a o\\m-s + \\uo\\h-s), 
\\ u \\l^(b^ s 2 ) - C ^ + e o){\\uo\\h-^ + TWaoWfji-s), 
Ml 1 ^- 5 ) + W u h 2 T (B^ 2 s ) ^ C {W u o\\h-s +T\\a \\ ti i-s). 
Here the constant C depends only on A, p,p,5,CQ. 

Proof. Set (3 = 1 — 6. Without loss of generality, we can assume that 771, 772, T small enough 
such that 

Cvi<\, C V2 <± CT< 1 -, (5.2) 

where C is the constant appearing in the following estimates and depending on X,p,,p, 6, cq. 
Step 1. Estimate for the transport equation 
Apply the operator Aj to the first equation of (|4,2p to obtain 

d t Aja + u ■ VAja = AjF + [u, Aj] • Va, 

with F = —(1 + a)divu. Making 1? energy estimate, we get 

1 d 1 /* 

- — ||Aja||l-- J 3 \A ja \ 2 divudx < (\\A j F\\ 2 + \\[u,Aj] ■ Va|| 2 ) || Aj-a|| 2 , 
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this gives 

rt 1 

\\A ja {t)\\ 2 < ||Ajao||2+ / (HAj-ii'lla + ||[«, A,-] - Va|| 2 + -||divu||oo||Aja|| 2 )dT, 

from which, we deduce that 



o 



I ^ 1 1 Too ( p>p \ < \\ao\\f,f> +C / ||divu|| _ 3 \\a\\~ Toa( ^ ,dr 

u t \ a 2,2> D 2,2 Jq ftP t \ a 2,2> 



J P,1 



+W F hl(Bl 2 ) + II^HKAj] • Va\\ L i {L2) \\ e2 . 



We get by Lemma 15.51 and 5 > 1 — | that 

ll-Fllrww 9 ^ < C\\ u \\ti ( rP+i\( 1 + 2 ). 

L T\ B 2,2> L T( B 2,2 ) V L^(B P l ) 

and by Lemma 15.71 

II^IIKA^-Voll^jll^ <C\\u\\ i +1 \\a\\ L¥{6 P y 

L T\ B p,l I 

Plugging them into (|5.3p and using (|5.ip . we infer from Gronwall's inequality that 



^ ) < C(||ao|| + H\tUBl+ 1 ))- 



L r (-^2,2 

Step 2. Estimate for the momentum equation 

Apply the operator Aj to the second equation of f|4.2|> to obtain 

d t Aju - div(TiVA j u) - V(Adiv Aju) = AjG + div([A i ,7i]Vn) + V([A, A^divu) 

and recall that 

G = -u-Vu - P P — Va + • Vu + 4-Vpdivn. 

p /r (r 

Making L 2 energy estimate(see the proof of Proposition 13. 3p . we obtain 

\\ u k ¥i B&) ^ ^(INII^ + IIGIIz^-!)) 

+C\\yP(\\[A j ,u]Vu\\ L i T{L2) + IKAj./llViillii^))!^, 

and 

-^TV- D 2,2 / T \ 2,2J 

+C||2^(T)(||[A„I7]V U ||^ (L2) + ||[A i5 7i]V U || LML2) )|| £2 . 
First of all, by Lemma |5.6( £ > 1 — -) and Lemma 12.81 we get 

||2^(||[A„!7]Vn||^ (L2) + ||[A„7i]Vn||^ (i2) )||, 2 



< C|| a IL 1 \\ u \\ti 

l t( B p,i) t ' 2 
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and 

\\2^ Wj (T){\\[^, V]Vu\\ Ll . m + \\[^Ti\^A^))\\p 

< C\\a\\„ 3 \\u\lf! /r/3+i v (5.8) 

Applying Lemma 15.51 and Lemma 12.81 to G, we obtain 

ll^llfi mf- 1 ) — C(II U IL a ||«||r2 f nj9 \ + a \\u\\fi ,f>i3+u 

Ll T\ t> 2,2 I v L^(B P 1 ) L T\ a 2,2l L^(B P ± ) T \ 2,2 / 

+ T II«HZ ?? K, 2 ))' ( 5 - 9 ) 

and 

+ T H«IIZ §? (^ 2 ))- ( 5 - 10 ) 

Here we need 1 — ~ < <j - 

Plugging (|5.7p and (|5.9|) into (|5.5p and using (|5.ip . we get 



II^Hl-^- 1 ) ^ Cdl^ll^ + T ll«lli ?? ( J B| i2 ) +^011^11^(51+1)), (5.11) 
and plugging (|5.8p and (|5.10p into (|5.6p . we have 

IN^j^t 1 ) + H u llz» (B£ a ) ^ C '(II«o|Ih-. + r||a|| £§?( ^ a) ) (5.12) 

Then Proposition 15. II follows by combining (|5.1ip - (|5.12p with (|5.4p . ■ 

5.2 Propagation of high regularity 

We will show that the solution constructed in Section 4 can also propagate the regularity of 
the initial data in Sobolev space with high regularity. 

Proposition 5.2 Assume that (p,u) is a solution of ([Qp on [0, T], which satisfies 

P>c , o€L ao (0,r;B| 1 ), u€L oo (0,r;s|7 1 )nL 1 (0,2 1 ;s|/ 1 ). (5.13) 

i/ (ao, uq) £ -fP x fP _1 for s > 3, then we have 

a € L°°(0,T;B1 2 ), u G L°°(0, T; S^ 1 ) n L^O, T; S^ 1 ). 

Proof. Due to (|5.13p . we can divide the time interval [0, T] into finite many small intervals 
[Tj, Tj + i] with i = 0, • • • ,N such that 

||a|| 3 < e, 

L-(T l ,T i+1 ;B p p 1 (^)) 

L^TuT^B^ ) L2(T l ,T i+1 ;B p p 1 ) 
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for some e small enough, see Remark 15. 31 Here the weight on each [Tj,Tj + i] is given by 

4 = J2 2k ~% 4 = (1 - e-^^-^ji 



^>fc 

First of all, we also have(see 15. 3|) 



IHIz~(Z| 2 ) ^ H a 0|U| 2 +C / H di HI f ll a llz~(Zf 2 ) dr 

+ II F IIZJ(^ 2 ) + ll 2JS H[ n ' A i] • Va|| L i (L2) ||, 2 . (5.14) 

By Lemma 15.81 we have 

ll^llzK^) - c ^ u hu&&) + IHL^fJMIz^ 1 ) + ^ a h r{ B h )\\ u \\ llAB i+i^ 

and by Lemma 15.101 

||2^||K A,] • Va|| z i (i2) ||, 2 < C(\\u\\ | +1 >llxr(BS, a ) + I^^^HI- . • f > 
from which and (I5.14p . we get by Gronwall's inequality that for £ 6 [0, Ti], 

IHIz~(%>) ^ C '(ll a oll^ 2 + ( 5 - 15 ) 
Making ZP _1 energy estimate for the momentum equation, we get (see [531 and (|5.6p ) 

'\ u k T{B &) < c'CIKoIU-i + ll^llzK^ 1 )) 

+C||2^(||[A,-,77]V«|| i i (i2) + lltA^/ZlVulLi^))!!,,, (5.16) 



and 

ll U IIZi(B 2 s + 1 ) + H U HZf(B|, 2 ) 

< ^(hob- + l|G|| Z x ( ^-x M) ) 

+C||2^(*)(ll[A,-,F]Vu|| itl(i2) + IICA^Vull^))^. (5.17) 

By Lemma 15.81 and Lemma 12.81 we have 



L t\ B 2,2 ( UJ )) \ L^(B P 1 ) L t\ B 2,2> Ll{B P 1 ) T\ B 2,2> 

+ fl L .3 It 7lfRS + l\ + ^ roo/6s \ "ft L — -1-1 "l ^ ft Too / R s 

L?°( B p,i^)) A 2 ' 2 ) * ( 2 ' 2) LU B p,i ) * ( 2 

and Lemma 15.91 and Lemma 12.81 

ll^'fjC*) (II [Aj, V«|| £ x (ia) + IKAj./ZlVttHrx^))!^ 



< C( a 3 llull j i/ns+h + I a| fooi-Ds Jul 3,, 
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which along with (|5.17p gives that 

W u ki(B&) + \\ u hi(B i2 ) < c{\\Mh*-i + ^kriB-iJML^i^ + *))> 

from which and (|5.15p . it follows that for any t E [0, Ti], 

H a IIZ-(B| 2 ) + MljiB^ 1 ) + \\ U k*(Bl 2 ) ^ C (W a o\\H° + IKHffs-i)- ( 5 - 18 ) 

On the other hand, by Lemma 15.81 and Lemma 12.81 again, we have 



H^llfiCR*- 1 '! < C Mlnll 3 ||it||f2/'R« \ + ||«|L a IMIfSfRs -i 

L t\ B 2,2 ) \ L'f(B P 1 ) t ( 2,2> L^(B P 1 ) L 'T < < tS 2,2> 



and by Lemma 15.91 and Lemma I2.8[ 



<C(|a| 3 ||ii|| ? i/rs+i\ + ||a|| ? tx)/Rs \\\u\\ 3,, 



which along with (15.16P and fj5. 15f) implies that for any t G [0,Ti], 

Then the proposition follows by repeating the above process in each interval [Tj,Tj+i]. I 

~ . 3 

Remark 5.3 Thanks to a G L°°(0, T; i?^ 1 ) ; one can choose JVi 6io enough so that 

2 ^ll A i«(*)lli-(0,T ; iP) < |- 

lil>iVi 

Fix iVj_, one can divide [0,T] inio finite many small intervals [Tj,Tj + i] mi/i i = 0, • • • , iV 
snc/i i/iai /or i = 0, • • • , iV 

£ 2 J lm|A i a(t)|| L ^ (TiiTi+i;L p ) < |. 
|j|<ATi 

T/ien we conclude that 

\\ a \\~ 3- 

L°°(T l ,T I+1 ; J B p p 1 (a;»)) 

<2 ^ 2 J 'f||A i a(t)|| LOO(T . iT . +1 . LP) + ^ 2 J 'fa;}||A j a(t)||L<»(T i ,T i+i; LP) < e- 
|i[>iVi |j|<JVi 
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5.3 Some technical lemmas 

i,i 1 

, J- ^ HI, q2,1 ^ w 

(a) if S2 < |, we have 



Lemma 5.4 Let p > 2,1 < qx,q2,q < oo with — + ^ = ^ . TTien i/zere ZioW 



f6j i/si < | — 1, we /lave 



W T f9\\~ s 1+S2 -i < 0\\f\\^ p i iiu)) \\g\\ z?( ^ 2 y 



(cj i/si + S2 > 0, we have 

llR{f,9) huB 3 2 \ +s ^hu>)) - cm HHK^)) M HH^ 2 y 

Proof. We first prove (a). Due to (|2.ip . we have 

Ai(T 5 /)= J] Aj(Sj'-igAff), 
b"'-il<4 

then we get by Lemma 12.11 and (|2.4ft that 

ll^/H +-a _ * = ||2^ (si+S2 -i ) o;,(r)||A i (r g /)|| i , (i2) ||, 2 

L A B 2,2 ( w )) 

< CII/llz^^H) 11911 ^ 2 ^)' 
where we used in the last inequality 



WSj-igW ^ <C Y ^||A, 9 || L , 2(L2) <C2^-^|b|| z , 2 ^ 2 



We next prove (b). Similarly, we have 

L t( B 2,2 M) 



< C 2 



J'(*i+aa-|) 



p> ujj (T) 1 1 Sj_ i / 1 1 L n (Lao ) 1 1 A jg 1 1 L « (i2) \\ i2 , 



and by Lemma 12,11 and ()2.4p , we have 

^o.7 \^ o^-l), . . I r r\ II a _ -f II . 

°) 



?3 



< C2^~^ Y 2 isl MT)\\A f J\\ L < THLP) 



< C2^-^\\f\U (Bai 
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which lead to (b). Finally, let us prove (c). Noticing that 

A J (R(f,g))= AjCAjv/A^), 

then we get by Lemma 12. II that 



+ s 2 - 



\W,g)\L q . n ... 

<C|| £ 2^^+^)a; J (T)||A J ,/|| L « (LP) ||A J , 5 || L?2(L2) ||, 2 
i'>i-3 

<c|| E E^^'^ii^'/ii^^jII^ii^^)!^ 
= II E E'"IU + II E E ■■■lU 

j'>j-3 e>j,j> i'>j-3 i'>£>i 

= 1 + 11. 

In the following, we denote {cj} by a sequence in £ 2 . Noting that 

2?~ e e e (T) < 2 j_j ' E 2J '~V r ) = 2 j ~ j 'wji{T), 
e>j,j' t>f 

we infer that 

I < C\\J2 ^Kr)2^ Sl+S2 )2^'||A J ./||^ 1(iP) ||A i ^|| i??(i2) 
i'>i-3 

< c\\ £ u; i ,(r)2(i-i')(^+^+i)2^ Cj ,||A / /|| L? x (iP) 
i'>i-3 



^2 " 3 Ltj, {iS 2 ,2> 



< C \\f\\l^{Bl^))\\9\\L^ 



and thanks to eg(T) < &j>(T) < toy(T) for I < j' , we get 

II < C\\ E ^"m2 j(si+S2) ||A,,/|| L , 1(iP) ||A J , 5 || i?2(i2) 

i'>i-3 

< C|| E "iK^Z^^^^IIA^/llpl ll^ll^^, 
i'>i-3 

^ C 1 1 / 1 1 1« (B^ («)) 1 1 9 1 1 L« 2 (B 2 S2 2 ) ' 

This proves (c). 

From Lemma 15.41 and its proof, it is easy to see that 

Lemma 5.5 Let p > 2, si < ~, S2 < | , Si + S2 > 0, and 1 < q, qx, q% < oo toit/j — 
T/ien we have 



\\fg\L 



s±+s 2 - 



< 



lub 22 n 



C \\fhl^B s \)y\\L q T 2 (B^ 2 



and if s± < | — 1, £/ien 



H/flL Q .s 1+ s 2 , 
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Lemma 5.6 Let p>2 and s € ( — |, |). Then there holds 

||2 J 1[A, J]V 5 || L , (L2) ||, 2 < CH/II a lbllzi ( ^)' 

L T\ B p,ll 

"T l°p,l( w )) 

Proof. The proof is almost the same as Lemma 12.91 We write 

[/, A.-jVg = [T/, A,-]V 5 + T A , Vi? / + R(f, A 3 Vg) - & s {T Vg f) - AjR(f, Vg). 

It follows from Lemma 15.41 (a), (c) that 

II^^IIA^Tv^lLi^H.^CII/ll s \\9\\ iHB s V) , 
p^^A^f^g)^^ <C||/|| f \\9\\l hi ^y 

Recalling (12. 5p and by Lemma 12.11 and (|2.4p , we get 

||^(T)2i s ||ri. V9 /|| L x (i2) ||, 2 

<C||a;,(T)2^^ 1 )||A i5 ||^ (L2) £ || A f f\\ L¥{Laa) || |2 

i'>i-2 

<C||2^ s+1 )||A j5 ||^ (L2) ^ ^(r)2^'l||A i ,/|| i?>(iP) ||^ 

i'>i-2 

< C||/|l 3 Hffllri f R S +i v 



Note that 



- Ty)dTd k h(2 :/ y)A j ,g(x - y)dy 



+2^ / h{V{x-y))d k S^ 1 f{y)A j ,g{y)dy, 
from which and Young's inequality, we infer that 

Ihcr^ifiy, A,]v 5 ||^ (i2) ||, 2 <c||/|| 3 

Summing up the above estimates yield the second inequality, the first one is similar. 
From the proof of Lemma 15.61 it is easy to see that 
Lemma 5.7 Let p > 2 and s G ( — |, | + 1). Then there holds 

||2^||[A J ,/]V 5 || L , (L2) ||, 2 < CH/II .l + A9k U B* 2 y 

t ^pa > 

The following lemmas are used in the proof of propagation of high regularity. Since 
proof is very similar to Lemma 15.51 -Lemma 15.71 we state them without proof. 
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Lemma 5.8 Let p > 2,s 1 ,s' 1 < + s% = + s' 2 > 0, and 1 < q,qi,q<2 < oo mi/i 

— + — = - . Then there holds 
qi 92 ? 



11/51 W 82 ") - c(ll/ll ^^) |bl1 ^ 2 ^) + 

and if si, s[ < | — 1, i/ien 

Lemma 5.9 Lei p > 2 and s > — |. T/ien i/iere ZioZds 

||2-||[A„/]V 5 |L, (i2) ||, 2 < ^dl/H^^^jllf H£5.(W + imiw l2 H)IML, (J £ } ). 

||2^(T)||[A,,/]V 9 || L , (L2) || £2 <C(||/L | jy^( W + ll/ll^>))ll^ fA f J- 
Lemma 5.10 Lei p> 2 and s > — | . T/ien i/iere /to/ds 

||2-||[A,,/]V,|| L , (L2) ||, 2 < ^(ll/IL^^^II^^) + 11/11^(^)11^1^^)- 

6 Hoff 's energy method 

In a series of papers (22J [231 EH GS], Hoff developed a method to construct global weak 
solution for the discontinuous initial data with small energy. An important property of this 
solution is that the density does not develop vacuum and can not concentrate if the initial 
density is bounded below and above. In this section, we will adapt his method to our case. 
Due to the stronger condition on the initial data, the restriction on the viscosity coefficient 
can be removed and the obtained estimate is also better. 
We set a(t) = min(l,t) and define 

A 1 (T)= sup ||Vn(t)||^ 2 + / f p\ii\ 2 dxdt, 
te\o,T\ Jo iR 3 



A 2 (T)= sup a(t) [ p\ii(t,x)\ 2 dx + / f cr|V 
tefo.Ti Jr. 3 Jo Jr. 3 



u\ 2 dxdt. 



= [0,T] 

Here and what follows, we denote 

f = ft + u-Vf. 

Theorem 6.1 Let (p,u) be a solution of satisfying 

p-peC([0,T};H 2 ), u£C([0,T};H 2 )nL 2 (0,T;H 3 ). 

There exists a constant £q depending only on A, fx, co,p such that if the initial data (po,uo) 
satisfies 

co < Po(x) < Cq 1 , x £ R 3 , 



/ \po(x) - p\ 2 dx + ||u ||m = C < e , 
Jr 3 
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then we have 



j < p(t,x) < 2CQ 1 , {t, x) G [0, T] x R 3 , 



Ax{T) + A 2 (T) <C*. 

Proof. Due to the assumption, there exists a < To < T such that the solution (p, u) satisfies 



co 



-? <p(i,x) < 2CQ 1 , (t,x) G [0,T ] x R 3 , 



i 



A^To) + A 2 (T ) <C 2 . 

Without loss of generality, we assume that To is a maximal time so that the above inequalities 
hold. In the following, we will give a refined estimate on [0, To] for the solution. We denote 
C by a constant depending only on A, fx, co,p. 
Step l.L 2 energy estimate 

[ \p(t,x)\u(t,x)\ 2 + G{p{t,x))dx + [ [ \Vu\ 2 dxdt 
Jr 3 2 Jo Jr 3 

< [ ~Po\u \ 2 + G( Po )dx < CC , (6.1) 
Jr 3 ^ 

where G(p) is the potential energy density defined by 
Due to (|1.4p . we have 

b(p)l < CG(p), (6.2) 

if 5 GC 2 with = </(p) = 0. 
Step 2. if 1 energy estimate 

Multiply (ll.ip by u and integrate the resulting equation over R 3 to obtain 

/ p\ii\ 2 dx = I (—u-X7P + pAu-u + \X7divu-u)dx. (6.3) 
Jr 3 Jr. 3 

By the continuity equation, we have 

d t P + div(uP) = divu(P - P'p). (6.4) 
Then integration by parts yields that 

/ -ii-VPdx= I (divu t (P-P(p))-(u-Vu)-VP)dx 
Jr 3 Jr 3 

= dt [ dwu(P - P(p))dx + [ ({P'p- P)(divu) 2 + Pd i u j d j u i )dx 
Jr 3 Jr 3 

<d t j divu(P - P(p)dx + C||Vii|||. (6.5) 
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Integration by parts again gives 



p / An • udx 
/r 3 



= — — dt I \Vu\ 2 dx -p I diU 3 di(u k d k u 3 )dx 
2 7r 3 Jr 3 

< ~% 9 t [ \Vu\ 2 dx + C [ \Vu\ 3 dx, (6.6) 



(6.7) 



and similarly, 

A / Vdivu • udx = —^-dt f \divu\ 2 dx — A f diU 3 d{(u k 'd k u 3 )dx 

JR 3 2 JR 3 JR 3 

< --dt [ \divu\ 2 dx + C [ \Vu\ 3 dx. 
2 Jb 3 Jr 3 

And by (JgTTJ) and ([63]) . we have 

/ divu(P - P(p)dx < ||divu|| 2 ||P - P(p)h < CC 1/2 1| divu|| 2 . 
7r 3 

Plugging (|63D-(I621) into ([631) yields that 

Ai(T ) < CC + ° y |Vu| 3 dxcft. (6.8) 



Step 3. H 2 energy estimate 

We take the material derivative to the second equation of (jl.lh to obtain 

pu t + pu-Vii + VP t + div(VP <g> u) 

= p[Au t + div(Aw ® u)] + A[Vdiv«i + div((Vdivn) <8) «)] 



Let us introduce some notations. We denote div(/ ® u) = Y2j=i dj(f u j)- For two matrices 

A = (oy)3 X 3 and i? = (foy)3 X 3j we use the notation ^4 : Z? = £)i 7=1 a y&ij an d is as usual 
the multiplication of matrix. 

Multiply the above equation by a(t)u and integrate by parts to get 

/ P\u\ 2 dx) - -a' I p\ii\ 2 dx 
V 2 Jr 3 y * 2 Jr 3 

= — I au J [djPt + 6xv(djPu)]dx + pa u- (Au t + div(Au ® u))dx 

JR 3 JR 3 

Vdivitt + div((Vdivu) g) ii) dx. 
Now integration by parts and use (I6.4j) to obtain 

aii 3 [djP t + dw(djPu)]dx = - [ a[(u-VP + P' pdwu)d jU j - d k u j djPu h ]dx 



+Xa / u 
/r 3 



'R 3 

= - / a[P' pdivudjii 3 - Pd k {u k dju 3 ) + Pd j (d k it i u k ) 

JR 3 

< Ccr||V«|| 2 ||Vn||2. 
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Integrate by parts again to obtain 



/ u ■ (Aut + div(Au ® u))dx = — u / [Vu : Vu t + u ® An : Vu] dx 
Jr 3 Jr 3 

Vu| 2 - V(u • Vu) : Vu + u (g) An : Vu 



< 



R j 



R ' 



R ' 



R 



|Vu| 2 - ((VnVn) + (n • V)Vn) : Vu - V(n • Vu) : V? 



|Vu| 2 - (VnVn) : Vu - div(Vu ® u) : Vu - (VuVu) : Vu - ((u • V)Vu) : Vu 

|Vu| 2 - (VuVu) : Vu + ((u • V)Vu) : Vu - (VuVu) : Vu - ((u • V)Vu) : Vu 
3/i 



|Vu| 2 + C|Vu| 



dx. 



To estimate the last term, note that 



div((Vdivu) (8) u) = V(u • Vdivu) — div(divuV <8) u) + V(divu) 2 , 
divu = divut + div(u • Vu) = divut + u • Vdivu + Vu : (Vu)'. 

Here A' means the transpose of matrix A. We have 



A / ii 
/r 3 

-A 
-A 



Vdivuj + div((Vdivu) (8) u) 



R 3 



R 



< 



R 



divudivut + divu(u • Vdivu) — divu(Vu)' : Vu + divu(divu) 2 

dx 

dx. 



d,x 



A 



|divu| 2 — divuVu : (Vu)' — divu(Vu)' : Vu + divu(divu) 2 
1 



-|divu| 2 + -|Vu| 2 + C|Vu| 4 



Summing up the above estimates, we obtain 

A 2 (T ) < CC + CA 1 (T ) + C f T ° [ a\Vi 

Jo Jk 3 

which along with (|6.8p gives 

/•To r rib r 

Ai(T ) + A 2 (T ) < CC + C / a\Vu\ 4 dxdt + C / \Vu\ 3 dxdt. 

Jo Jn 3 Jo Jn 3 



7 u\ dxdt, 



(■T {) 



(6.9) 



(6.10) 



It remains to bound the right hand side. For this purpose, we set 
F = + X)divu- P(p) + P(p), oj = Vxu. 

By (I6.4j) . we have 

d t (p-p(p)) + u-v(p-p(p)) 
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which implies that 

c\\p{ P ) - p{p)\\i 

<d t J{p- p(p)fdx + s\\p( P ) - p(p)\\i + c(\\f\\1 + l|v«||2 

Taking 5 small and integrating the above equation on [0,To], we get 

rT f a(T) 



/ a\\P(p) - P(p)\\Ut < C sup \\P(p) - P(p)\\ls + C / \\P(p) - P(p)\\ 

JO t€[0,T] JO 

/•To 

+C a\\F\\idt + CC 
Jo 

< CC + C [ ° a\\F\\ 4 L4 dt, 
Jo 



~^%dt 



from which and Lemma 16.21 we infer that 

I a\Vu\ 4 dxdt 
o Jr. 3 

T rTo 



<C I ° o{\\F\\i + \\u>\\\)dt + C / ° a\\P(p) - P(-p)\\\dt 
Jo Jo 

<C t ° a{\\Vu\\ 2 + \\P(p) - P(p)\\ 2 )\\fmf 2 dt + CC 



o 



and 



<C sup (a 1/2 ||p 1/2 n|| 2 (||V U || 2 + C 1/2 )) \\p 1/2 u\\ 2 2 dt + CC 
te[o,T ] Jo 

< C(A{ /2 + C 1/2 )4 /2 ^i(T) + C < CC , 



To r 

/ \Vu\ 3 dxdt 
JR 3 

r&(To) r rTo r 

< / \Vu\ 3 dxdt + / / \Vu\ 3 dxdt 

JO Jr 3 Jo-{T ) JR 3 

&(To) 3 3 3 f T 



<[ l|V«|||(||/ou||| + \\P - P(p)\\l)dxdt + [ ° / (|Vu| 4 + \Vu\ 2 )dxdt 
JO Ja(T ) Jr. 3 

3 r°(To) 3 

<C sup \\VuM \\pu\\l+CC 

te\0,a(To)\ JO 



< CA\ + CC < CC 4 . 
Plugging them into ()6.10j) and taking e small enough depending on p, \,p~,co, we conclude 

Atpo) + A 2 (T ) < (6.11) 
Step 4. Lower and upper bound of the density 

26 



Set L = log(p), which satisfies 

(2p + X)L+(P(p)-P(p)) = -F 

For < t < a (To), we have 

/ \\F\\ L ~dt < C / \\F\\l 6 \\VF\\l 6 dt 
Jo Jo 

< c r To \\pu\\i 2 \\ P u\\i 6 dt 

Jo 

f<r(To) i i 

< c / Wp^m.wvuWl.dt 

Jo 



JO 

< ccl 

which implies that for t < a (To), 

inf(logpo(^)) - CC\ -Ct< log p(t,x) < sup(logp (a:)) + Cc| + Ct, 
hence, one can choose £q,t sufficiently small such that for t < r < cx(To), 

^c <p(t,x) < j^o 1 . (6.12) 

For t < t\ < t% < To, we have 

2 ||F|| L «,cft < C / 2 Hp^tilltllVtillitft 



L 2|| v u,|| L2 



< C(t 2 -t x fl 2 A||p 1/2 u||| 2 + ||Vu||| 2 )* 

Jti 

< »7(* 2 - *i) + CC 

for any r\ > 0. Now let T* be a maximal time so that (|6.12p holds. If T* > T, then we are 
done. Otherwise, we have 

p(T\x(T*)) = \c Q l or p(T\x(T*)) = \c , 
where the curve x(t) is defined by 

— x(t) = u(t,x(t)), x(0) = x. 
If p(T* ,x(T*)) = fcQ 1 , then there exits T* > T\ > r such that 

p(T u x(T x )) = Cq 1 , Cq 1 < p(t,x(t)) < ^c 1 for t G [T^T*]. 
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Then we conclude by (|1.4p that 

L(T* ,x(T*)) — L(T\,x{Ti)) < CCo, 
if we choose r\ small enough. If we choose £o sufficiently small, this in turn implies that 

p{T\x(T*))<\c Q \ 
which contradicts with the definition of T*. Similar argument yields that 

p(T*,x(T*)) >^c , 
if p(T*,x{T*)) = |c . Hence, we have 

^c <p(t,x) < -c^ 1 , for 0<i<T . (6.13) 
With (|6.1ip and f|6. 13j) . we conclude the proof by continuity argument. ■ 
Lemma 6.2 For any p 6 [2,6], there exists a constant C depending only on p,X such that 

\\Vu\\ p < C{\\F\\ P + \\lj\\ p + \\P(p) - P(p)\\ P ), 

5=£ Sp-6 

\\Vu\\ p < C\\Vu\\ 2 2p (\\pu\\ 2 + \\P(p)-P(p)\\ 6 ) 2 p . 
Proof. By the elliptic estimate and the definition of F, we have 

llVullp < C(||divn|| p + \\lo\\ p ) < C(\\F\\ P + \\co\\ p + \\P(p) - P(p)\\ P )- 
Noting that 

AF = div(pu), pAoo = V x (pii), 
we have by the elliptic estimate that 

||VF|| 2 + ||Vw|| 2 < C\\pu\\ 2 . 
Then by Holder inequality and Sobolev inequality, we get 

6— p 3p— 6 

\\Vu\\ p < C||Vu|| 2 2p ||Vw|| 6 2p 



6—p 



3p-(> 



< C\\Vu\\ 2 ^ (\\F\\ 6 + ||u;|| 6 + \\P(p) - P®h) 2 * 



6—p 



3p-6 



< C\\Vu\\ 2 2 * (\\pu\\ 2 + \\P(p)-P(p)\\ 6 )^ 
This finishes the proof of lemma 
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7 Continuation criterion 

In order to extend a local solution to a global one, we need to establish a continuation criterion 
of smooth solution. The proof is motivated by [33]. 

Theorem 7.1 Let (p,u) be a solution of U.l\) satisfying 

p(0)>0, p-p€ C([0,T};H 2 ), u G C([0, T];H 2 ) n L 2 (0, T; H 3 ). 

Let T* be the maximal existence time of the solution. IfT* < +oo, then it is necessary 

limsup (||p(t)||oo + |K*)|| g ) = +oo (7.1) 
tyr* 

for any q > 3. 

Remark 7.2 Compared with Theorem 1.3 in [33], the restriction on the viscosity coefficient 
(A,/i) is removed, but the price to pay is to impose the condition on the velocity u. Due to 
p(0) > 0, the compatibility condition on the initial data is also removed. 

Proof. We use the contradiction argument. Hence, assume that T* < oo and 

sup (||p(i)||oo + K*)|| 3 ) =M<+oo. (7.2) 
te[o,T*) 

In what follows, we denote C by a constant depending on T,M, \\uo\\ H 2 and ||po ~7>\\h 2 - 
First of all, we have 

I -p{t,x)\u(t,x)\ 2 + G{p{t,x))dx + [ [ \Vu\ 2 dxdt<C. 
Jr. 3 2 J j R 3 

This gives by (j7.2j) and (j6.2j) that for any r G [2, oo], 

\\P ~ PlU°°(0,T;£'-) + \\VP~u\\l™(0,T;L 2 ) + II ^u\\l 2 {0,T;L 2 ) < C. (7.3) 

Let v = L^ 1 VP(p) be a solution of the following elliptic system 

Lv = pAv + XVdivv = VP(p). 

By elliptic estimate and (|7.3p . for r £ [2, oo) 

||Vu|| r < C\\P(p) - P(p)\\r < C, \\V 2 v\\ r < C\\Vp\\ r . (7.4) 

Introduce a new unknown w = u — v. Note that divu> = (A + p)divu — P(p) + P(~p) is a 
so called the effective viscous flux introduced by Hoff [22J. It is easy to find that w satisfies 

pdtw — pAw — AVdivw = pF, 

with 

F = -u ■ Vu + L' 1 Vdiv(p(p)u) +L' 1 V{(pp'{p) -p(p))divu). 
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Multiplying the above equation by dtw and integrating by parts, we obtain 

-/ /x|V^| 2 + A|divu;| 2 a!x + - / p\d t w\ 2 dx < -\\y/pF\\l 
at Jr3 2 J R 3 2 

By Sobolev inequality, interpolation inequality, (|7.2p and (|7.4p . we get 

Hy^u • Vit||i2(Q) < C||tt||, 7 ||V , u|| 2 g < C(||Vu>|| 2 9 + [|V«[| i q ) 

q — 2 q — 2 q — 2 

< C[|Vtw[| 2 + e||V 2 u>|| 2 + C 

< C + C||Vu|| 2 + e||V 2 uj|| 2 

for any e > 0, and by elliptic estimate, (|7.3p and (|7,4p . 

\\^pL- l Vdiv[p(p)u]\\ 2 < C\\p(p)u\\ 2 < C\\y/pu\\ 2 < C, 
and by Sobolev inequality, 

W^fpL^V (pp' -p)div«|| 2 

< II/ ~ P||3||L _1 V(pp' - p)divw|| 6 + C\\(p - p)divu|| e + 
<C||V«|| 2 + C*||^|| 2 <C + C||Vn|| 2j 



here we used 



Hence, 



p 2 / \u\ 2 dx<2 \p-p\ 2 \u\ 2 dx + 2 p 2 \u\ 2 dx. 
Jr 3 Jr 3 Jr 3 



(7 



II VpF\\ 2 < C(l + [|V«|| 2 ) + e||V 2 u>|| 2 . 
Noting that Lw = pdtw — pF, we get by elliptic estimate that 

||V 2 u;|| 2 < C(\\pd t w\\ 2 + \\pF\\ 2 ) < C(\\^-pd t w\\ 2 + \\yfpF\\ 2 ), 
which implies by taking e small enough that 

llv^lh < c(i + l|v«y + \\\y/pd t w\\ 2 . 

This along with (|7.5p gives 

l|Vui|| LO o (0iT;L 2) + IIV^^^IIl^o.t^ 2 ) + llv 2 ^!!^^.^) < c, 

which in turn implies by (|7.4p and Sobolev inequality that 

l|V«|Uoc (0 ,r;^) + UVull^p^ < C for rG [2,6]. (7 

Now let us turn to the high order energy estimate. Thanks to the proof of Step 3 
Theorem 16.11 we have 



4" / tp\u\ 2 dx + pt [ \Vu\ 2 dx + Xt [ Idivul 

at j R 3 y R 3 y R s 



2 dx 



<C + Ct \Vurdx + C / \Vu\ 6 dx. (7- 
/r 3 Jr 3 
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Noting that 

fiAw + AVdivii; = pit, 

elliptic estimate yields that 

||V 2 W || 2 < C||Hl2 < C||v^u||2. 
Then by Sobolev inequality, (|7.4|) and (|7.6p . we infer that 

||V«||| < ||V«|| 2 ||Vu||| < C||V«|||(||Vto||e + ||Vu|| 6 ) 

< C7[|V«||1(1 + ||V 2 u;|| 2 ) < C||V«||1(1 + W^puh), 

3 3 3 

||Vu||i < ||Vu||f ||Vu||£ < C7||V«|||. 

Plugging them into (|7.7p and noting ||Vu(i)||6 G L 2 (0,T) by (|7.6p . we deduce by Gronwalbs 
inequality that 

fT 

|2,v„, 1 / / j.|v7„\ |2, 



/o|i»rda;+ / / t|Vitrdxdt < C. (7.8) 
r 3 j Jn 3 

from which, elliptic estimate and Sobolev inequality, it follows that 

l|V 2 H|z,2( 0!T;Z ,r-) < C\\pu\\ L 2( 0jT . Lr) < C for r€ [2,6]. (7.9) 

To obtain the high order estimate for u, we need to use the continuity equation. Take the 
derivative with respect to x for the first equation of (jl.ip to obtain 

d t Vp + {u ■ V)Vp + VuVp + divuVp + pVdivu = 0. 

For r £ (3, 6], making L r energy estimate yields that 

\Vp\ r dx 

C I \Vu\\Vp\ r dx + C f plVdivulNpl^dx 

JR 3 Jr 3 



d_ 

dt 



R 3 

< 

Jn 3 ' " Jr 3 

^CllVnlUllVpll^ + CHV^II.IlVplir 1 

2„..ll 1 iiv72„,ii MlT7„lir-l 



< CdlVwIU + ||V^|U)||Vp||; + C(||V 2 u;|| r + ||V 2 r,|| r )||Vp|| 
Then by (|7.4p and Sobolev inequality, we get 

^ / \V P \ r dx < C{1 + \\VvU + \\V 2 w\\r) ||Vp||;, 
at j R 3 

and by Lemma 17.31 and elliptic estimate, 

||Vu||oo < C(||Vu|| 2 + \\Vv\\^ Loa \n{e + \\Vv\\ w i, T )) 
< C(||p-p|| 2 + ||p|Uhi(e+||Vp|| r )), 
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from which and (|7.9p . we deduce by Gronwall's inequality that 

||Vp|| r < C, (7.10) 

which along with (|7.9p and (|7.4|) gives 

II V2 ^IIl 2 (0,T;L'-) < l|V 2 tf || L 2 (0iT;i r) + ||V 2 u|| L 2( 0)T . L r-) < C. (7.11) 

This in turn implies (|7.10p also holds for r € [2, 6]. 

Now we are in position to give H 2 estimate of (p, u) . First of all, we have 

j t \\p - n 2 m < c ( l + l|Vn|U)||p - p\\m + C||V 3 «|||. (7-12) 

While by (|7JD|> . 

\Mfm)\\ 2 < ||Vp|| 3 |H| 6 + C||V«|| 2 < C||V«|| 2 , 
this along with elliptic estimate and (|7.8p implies that 



Indeed, we have 



N 2 w(t)\\t+ t \\V 3 w(T)\\tdT < C. 
Jo 



' \\V 3 w(r)\\ 2 2 dT < \\V 3 w(r)\\ 2 dr + f \\V 3 w(r)\\ 2 dr 
J o Js 

t 

£ C + 5- 1 I tW 3 w(t)\\Ut < C. 





Hence by (|7.4|) . we get 

IIVVOH 2 . + / \\V 3 u(T)\\ 2 2 dT <C + C f ||V 2 p(r)|||dr. (7.13) 
Jo Jo 

Summing up ()7.12j) and (|7.13p . we conclude by Gronwall's inequality that for < t < T*, 

\\u(t)\\ H Z + \\ P (t)-p\\ H 2 <C. 

This ensures that the solution can be continued after t = T*. ■ 
Lemma 7.3 Let r > 3 and u £ W 1,r . Then there holds 
IMloo < C\\u\\2 + C\\u\\ a ln(e + 

00 ,00 

Proof. The proof is standard. For the reader's convenience, we present a proof. We decom- 
pose u into 

N 

u = S-nu + Aju + A? M ' 

j=—N j>N+l 
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from which and Lemma |2. 11 we infer that 

Hoc < 2-l N \\u\\2 + (2N + l)\\u\\zo +Cy"2^'||A 



j>N 

< 2-l N \\u\\ 2 + (2N + l)\\u\U +C2^ 1 -^ N \\u\\ m , r . 

oo ,00 

Taking N £ N such that 

2^ 1 -^ N \\u\\ wl ,r ~ 1, 

the desired estimate follows easily. ■ 

8 Proof of Theorem 11.11 

Theorem 14. 1\ Proposition 15.11 and Proposition 15. 21 ensure that there exists a solution (p,u) of 
(jl.ip satisfying 



co 



< p < 2cq\ p-pe L°°(0, T; n B| >2 ), 



« e i°°(o, r ; * n s^ 1 ) n ^(o, r ; b; x n s^ 1 ). 

Moreover, there hold 

II^-^IIl-^ 1 -*) ^ C'CllPo - pIIhi-s + ||«o||jj-«)i (8-1) 

ll^llz-cs-l) <C(l + E )(\\u \\^ s +T\\p -p\\H^), (8.2) 

H n llz^(B 2 2 -/) + IMIz^bJ/) ^ CCIKII^ +T||p -pII^i-*), (8.3) 

INI 3 , , + llull s < r/ 2 < 1. (8.4) 

4(B* + ) L^B^) 

and by Remark 15.31 the existence time T has a lower bound 

T ~ c(i + \\ P0 -p\\.i )<"+*>/(-§> =Tl ' 

Here and what follows, the constant C depends only on A, p,co,p,p, s. Due to ||uo||h-« < C2, 
we get by ([12]) and ([83]) that 



ll u lli°°(B- 5 ) - C(l + E ) (c 2 + T\\p Q - p\\fji-s) 

\\u\\?i /o2-«s + ||u||r 
• £v tI-"2,2 J 

Then for r £ (1,2), we get by Lemma 12.11 that 



< cdhllz-cB^) + IMIz^*)) 

<C(1 + £ )(C2 +THP0-PIIH1-*), 
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and similarly, 

IMIz4(L2) < c^lhllr-^ + IIVniiir^)) 

< C{1 + Eo){c 2 + T\\po-p\\h^)- 
Hence, there exits a point to £ (0, T) such that 

||«(to)||ffi <C(l + ^,)(c2/T§+Ts ll/oo -p||^_ 4 ). (8.5) 
Recall that the density p satisfies 

<h(p — p) + u- V(p — p) = -divup. 
Making L 2 energy estimate, we get 

dt\\p-ph < l|divu||oo||/9 - p|h + C||divu|| 2 . 
Then by Gronwall's inequality and (|8.4|) . we get 

||/o(t)-p|| 2 < C(||po-p||2 + T5||Vn|| i2T(i2) ) 

< C(ci + T3||uo|| 2 ). (8.6) 

Here we used 

II v, "IIl2(l2) < C(IIPo -ph + IKIh), 
which follows from the energy inequality (|6.ip . Set 



.2 



^6 



C(l + \\po - p\\ a +||uo||2) 
and T* = min(Ti,T2). We choose ci,C2 small enough so that 

Cci<^>, c(l + E )c 2 /rf 

i.e, c 2 < e r|/(4C(l + E )). Then it follows from JSSK and (USD that 

||«(*o)||2 + [|it(to)[|fri < £<)• 

Theorem 16.11 ensures that 

j<P<4c -\ \\u(t)\\ L¥{Lq) <C. 
So, the solution can be extended to a global one by Theorem 17.11 
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